Mon. Not. R. Astron. Soc. 000, pHlS^ (2001) Printed 1 February 2008 (MN WT^ style file vl.4) 



Application of a relativistic accretion disc model to X-ray 
spectra of LMC X-1 and GRO J 1655-40 



o 
o 



(N 
(N 



Marek Gierlinski^'^, Andrzej Maciolek-Niedzwiecki^ and Ken Ebisawa^ 

^Astronomical Observatory, Jagiellonian University, Orla 171, 30-244 Cracow, Poland 
'^N. Copernicus Astronomical Center, Bartycka 18, 00-716 Warsaw, Poland 
'^Lodz University, Department of Physics, Pomorska 149/153, 90-236 Lodz, Poland 

■^Code 660.2, Laboratory for High Energy Astrophysics, NASA/Goddard Space Flight Center, Greenbelt, MD 20771, USA 
(also at Universities Space Research Association) 



Accepted on 2001 March 21 



ABSTRACT 



> 

(N 

m 
m 
o 



o 



We present a general relativistic accretion disc model and its application to the 
soft-state X-ray spectra of black hole binaries. The model assumes a flat, optically thick 
disc around a rotating Kerr black hole. The disc locally radiates away the dissipated 
energy as a blackbody. Special and general relativistic effects influencing photons 
emitted by the disc are taken into account. The emerging spectrum, as seen by a 
distant observer, is parametrized by the black hole mass and spin, the accretion rate, 
the disc inclination angle and the inner disc radius. 

We fit the ASCA soft state X-ray spectra of LMC X-1 and GRO J1655-40 by this 
model. We find that having additional limits on the black hole mass and inclination 
angle from optical/UV observations, we can constrain the black hole spin from X- 
ray data. In LMC X-1 the constrain is weak, we can only rule out the maximally 
rotating black hole. In GRO J1655-40 we can limit the spin much better, and we find 
0.68 < a < 0.88. Accretion discs in both sources are radiation pressure dominated. 
We don't find Compton reflection features in the spectra of any of these objects. 
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1 INTRODUCTION 

Black hole binaries (BHB) are generally observed in one of 
the two spectral states, as determined by their X-ray/7-ray 
energy spectra. In the hard state most of energy is radiated 
above ~ 10 keV. The spectrum is dominated by a power 
law with the photon spectral index F < 2.0 and a high- 
energy cutoff above 100 keV (Grove et al. 1998). A common 
interpretation of this emission is in terms of Comptoniza- 
tion of soft seed photons by thermal, hot, optically-thin elec- 
tron plasma. Cyg X-1 (Gierlinski et al. 1997) and GX339-4 
(Zdziarski et al. 1998) are typical examples of the hard state 
BHB. In the soft state most of the energy is radiated in the 
soft X-rays, below ~ 10 keV. This energy range is dominated 
by a blackbody-like component with characteristic temper- 
ature ~ 1 keV. This component is usually attributed to the 
thermal emission of a cold, optically-thick accretion disc, 
extending down to the marginally stable orbit (Shakura & 
Sunyaev 1973) . The disc spectrum is often accompanied by a 
high-energy tail, which can be described as a power law with 



a typical photon spectral index F ~ 2.0-2.5, extending into 
7-rays without apparent break or cutoff (Grove et al. 1998). 
LMC X-1 (Schlegel et al. 1994), LMC X-3 (Ebisawa et al. 

1993) , GS 1124-68 just after the outburst (Ebisawa et al. 

1994) or Cyg X-1 between May and August 1996 (Gierlinski 
et al. 1999) are the examples of the soft state BHB. 

At first approximation the soft component in the soft 
state can be described by a multicolour disc model (here- 
after MCD; Mitsuda et al. 1984), which spectrum is a 
sum of blackbodies with radial temperature distribution 
r(r) (X r""^/*. This model has been commonly used for 
spectral fitting by many authors. Though MCD model ap- 
proximates the disc spectral shape well, it ignores the in- 
ner torque-free boundary condition and parameters derived 
through it are incorrect. An improvement of the MCD model 
is a multicolour disc model in the pseudo-Newtonian poten- 
tial (Gierlinski et al. 1999) , taking properly into account the 
boundary condition. 

In vicinity of the black hole relativistic effects become 
important. The effects of relativity on the emerging disc 
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spectrum have been studied by Cunningham (1975), Laor, 
Netzer & Piran (1990), Asaoka (1989), Hanawa (1989), Fu 
& Taam (1990), Yamada et al. (1994) and others. Hanawa 
(1989) calculated disc spectra around a non-rotating black 
hole, creating a fast-computing model that have been used 
for spectral fitting several times (e.g. Ebisawa, Mitsuda & 
Hanawa 1991; Makishima et al. 2000). However, this model 
is limited to Schwarzschild metric only and effects of light 
bending have been neglected. More advanced model was de- 
veloped by Laor et al. (1990), who studied accretion discs 
around rotating black hole, taking into account both rel- 
ativistic effects affecting radiation and vertical structure of 
the disc. Laor (1990) applied this model to the AGN spectra. 

In this paper we show an application of the general rel- 
ativistic (hereafter GR) disc model to the soft state X-ray 
spectra of BHB. We focus on estimating the black hole spin. 
The model thoroughly treats relativistic effects in the deep 
gravitational potential of the black hole. In Section ^ we give 
description of the model and show examples of the model 
spectra. Then, we fit the X-ray spectral data of LMC X-1 
(Section |^) and GRO J1655-40 (Section |!|) with the GR 
model. In Section H we discuss the obtained results. 



2 MODEL DESCRIPTION 

We consider an accretion disc around a black hole (e.g. Page 
& Thorne 1974). The model is based upon certain assump- 
tions. We assume that the disc rotates in the equatorial 
plane around the Kerr black hole; that it is in a steady 
state; that it is axially symmetric; that it is fiat and geomet- 
rically thin; that it is optically thick; that it locally radiates 
away the dissipated gravitational energy as a blackbody; 
that there is no emission below the marginally stable orbit. 

A black hole is characterized by its mass, M, and angu- 
lar momentum, J. We use dimensionless spin a = Jc/GM^ , 
where < a < 1 and a — corresponds to the non-rotating 
(Schwarzschild) black hole. In accreting systems the radi- 
ation emitted by the disc produces a counteracting torque 
and the black hole cannot be spun-up beyond a — 0.998 
(Thorne 1974). Therefore we do not exceed this value in our 
model and refer to it as to the maximally rotating black 
hole. The efficiency of accretion, rj, varies from -q ~ 0.057 for 
the non-rotating black hole to « 0.32 for a = 0.998. 

The spectrum of the disc, _Fy, seen at infinity is com- 
puted by means of a photon transfer function, T, which 
describes travel of photons from the point of origin to the 
distant observer. The photon of frequency is emitted at 
radius r^, at cosine angle fie = cos9e and then perceived at 
frequency Vo by the observer at cosine angle fio = cosi. The 
observed fiux is 



'i'o(Mo) = ("^) J '^dcff J dre J d^e 

X rT{a,re,fJ.e,gcff,fJ,o)N{re,fJ.e, — )• 

9cS 



(1) 



Subscript 'e' denotes quantities measured in the local frame 
CO- moving with the disc and subscript 'o' denotes quanti- 
ties measured by the observer at infinity. N(re, fJ,e, = 
Nv^ife, fie) is the specific photon number intensity of the 
disc emission, geff = Vo/ve is the effective redshift of the 
photon, D is the distance to the source and Rg = GM/c? is 



the gravitational radius. Numerical constants are included 
in the transfer function. 

The transfer function treats the special and general 
relativistic effects affecting the spectrum. It takes into ac- 
count the Doppler energy shift from the fast moving mat- 
ter in the disc, the gravitational shift and the light bend- 
ing near the massive object. It includes integration over the 
azimuthal angle (j>e. Each element of the transfer function, 
T{a,re, fJ.e, gcS, fJ,o), was computed by summing all photon 
trajectories for all angles ijie, at given (re, fie), for which re- 
quired (goff, Mo) were obtained. The details of the transfer 
function computation are given in Appendix ^ A similar 
approach to computing the black hole disc spectra has been 
applied, e.g., by Laor et al. (1990) and Asaoka (1989). Note 
that our construction of T is slightly different than that of 
Cunningham (1975), as the latter one involves the relation 
of the specific intensities at the emitter and the observer, 
respectively, given by Liouville's theorem, while our model 
is based on counting individual photons. 

The local gravitational energy release per unit area of 
the disc, per unit time is (Page & Thorne 1974) 



Qix) = 



SMjC^ 1 

SttG^M^ a;4(a;3 ~ 3x + 2a) 

3(xi — a)^ 



X — Xo 



2 Va;o/ 



Xi{xi - X2)(xi — 3:3) 

3(2:2 — a)^ 

X2{X2 - Xi){x2 - X3) 

3{x3 — a)^ 

X3{X3 — Xl){x3 — X2) 



X — Xi 



xo 


— Xl 


X 


-X2 


Xo 


— X2, 


X 


- X3 



, xo — X3 , 



(2) 



where x = r^"^ , xo = r]l^ , xi = 2 cos(i arccos a — 7r/3), 
X2 ~ 2 cos(i arccos a -I- 7r/3) and X3 — —2 cos(| arccos a). 
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is the marginally stable orbit radius, where yli = 1 + (1 



a 



2x1/3 



[(l-fa)i/3 + (i_a)i/3] and A2 = {3a^+Ai)^/^. Both 
re and Vms are expressed in units of Rg. 

The photon number intensity, N,^^, can be derived 
from the energy release rate, Q. However, the relation 
between these two quantities is not simply Q{re) = 
J hveNv^{re, fie)AveAfLe, as assumcd in previous similar cal- 
culations (e.g., Laor et al. 1990). The photon number in- 
tensity which is used in equation (|^), is defined in terms of 
the distant observer coordinate frame. On the other hand, 
Q is defined as the energy dissipated per unit proper time 
per unit proper area, as measured by the observer orbiting 
with the disc. The corresponding dissipation rate measured 
by the distant observer will be affected by kinematic and 
gravitational time dilation and length contraction. We em- 
phasize that these effects are not included in the calculation 
of the single photon trajectory, and are treated separately 
when applying the transfer function formalism. 

In Appendix |b| we derive transformation of the time 
and the disc surface area between the disc rest frame and 
the frame of the distant observer: 



At' = PtAt, 



Ar' A<j)' — PsArA(j), 



(4) 



where {t, r, cf>) and {t' , r' , (j)') are the coordinates of the refer- 
ence frame of the distant observer and the disc, respectively. 
In order to find a locally emitted spectrum we would 
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need to compute the vertical structure of the disc. However, 
this problem did not find a satisfactory general solution yet. 
Instead, we simply assume that each point of the disc radi- 
ates like a blackbody and introduce two corrections. First, 
we note that the Thomson scattering dominates over ab- 
sorption in the inner part of the disc and the local spectrum 
is affected by Comptonization. At a given radius, Ve, we ap- 
proximate the spectrum by a diluted blackbody. 



io-« 



-S'^ = -F^B^^(fcoiTcft), 

/col 



(5) 



where i?^^ is the Planck function, /coi = Tcoi/Tcff is the spec- 
tral hardening factor, Tcoi is locally observed colour temper- 
ature and TcEF is the effective temperature, related to the lo- 
cally dissipated energy, Q, as uT^g — Q. Following Shimura 
& Takahara (1995) we use /coi = 1.7. 

The second correction concerns angular distribution of 
local emission. We assume a linear limb darkening in form: 



r 1 + <5Me „db 



(6) 



where S = 2.06 in classical electron scattering limit (e.g. 
Phillip & Meszaros 1986). 1^^ is the specific intensity in 
the local frame, co-rotating with the disc. Factor Iq is 
found from the requirement of energy conservation. Namely, 
the power emitted by a limb-darkened surface element dA, 
2n Iiy^ (/ie)Med/iedA, should be equal to the power emit- 
ted by the same surface element without limb darkening, 
■kB^^ dA. From this we find 



3 + 25' 



and 

If, (Me 



1-1- S^e „db 
■3 „ „^ B,. . 



(7) 



(8) 



3 + 26 

Then, the specific photon number intensity of the disc, 
measured by the distant observer, is given by 



(9) 



There is, however, one effect that we have not taken into 
account. Due to gravitational focusing some of the emitted 
photons return to the disc, where they are reprocessed or 
scattered. Since the returning photons alter the locally emit- 
ted spectrum recurrently, taking them into account would 
require enormous computing time and would make spectral 
fitting practically impossible. In our model returning pho- 
tons are lost, so the disc luminosity and temperature are 
underestimated for the spin close to maximum. This effect 
was studied by Cunningham (1976). He found that return- 
ing radiation is negligible for a < 0.9. The difference between 
the actual locally generated flux and that in the absence of 
returning radiation is only a few per cent over most of the 
inner disc for the maximally rotating black hole. At lower 
spins, the effect is much less pronounced. In this paper we 
compute the GR disc spectra for a = 0, 0.25, 0.5, 0.75 and 
0.998. Though the a = 0.998 model underestimates temper- 
ature and luminosity, the interpolation between the models 
would give accurate results for a < 0.9, and slightly overes- 
timated spin measurement otherwise. 

Figure |l] shows the GR disc spectra for different val- 
ues of the black hole spin. We clearly see how the accretion 
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Figure 1. Spectra from the GR disc model for a = 0, 0.25, 0.5, 
0.75 and 0.998 (from the bottom to the top). The disc extends 
from the marginally stable orbit to Tout = lO**. The spectra were 
computed for a IOMq black hole at 1 kpc, with accretion rate 
A/d = 10^* g s"-'^, disc inclination angle i = 45°, /coi = 1.7 and S 
= 2.06. 
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pseudo-Newtonian disc 

relativistic disc, 5 = 
relativistic disc, 6 = S.06 




E (keV) 

Figure 2. Comparison of the GR and PN disc models for differ- 
ent inclination angles. For both discs r^^ = 6. The spectra were 
computed for a non-rotating IOMq black hole at 1 kpc, with ac- 
cretion rate = 10^* g s~^, /coi = 1.7. The effect of the limb 
darkening in the relativistic disc is also shown. 
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efficiency grows witti increasing spin, when energy is dissi- 
pated deeper in the gravitationai potential. On Figure ^ we 
compare spectra of the GR modei (with a = 0) and the 
pseudo-Newtonian (PN) modei (Gieriihslci et ai. 1999). The 
PN model approximates temperature distribution with good 
accuracy, but effects of Doppler and gravitational shifts, 
light bending and focusing are neglected and light propa- 
gates through flat space. With increasing inclination angle 
both energy of the peak and normalization of the relativistic 
disc spectrum (dashed curve) increase in compare to the PN 
disc spectrum (dotted curve), mostly due to Doppler shifts 
and increase of the apparent disc area. A similar effect was 
found by Fu & Taam (1990). On the same figure we also show 
the effect of the hmb darkening (solid curve) , enhancing the 
spectrum for lower and diminishing it for higher inclination 
angles. Next, we check how the parameters obtained with the 
PN model relate to these from the GR model. To do this, 
we created fake ASCA spectra using the GR model with the 
following parameters: M = lOM©, a = 0, Md = 10^* g s"\ 
D = 1 kpc, Tin = 6, Tout — 10*, /col = 1.7. The spectra were 
computed for (5 = and 2.06, for several inclination angles. 
Then, we fitted the PN disc model to the spectra finding 
the black hole mass and the accretion rate. The results are 
presented on Figure ^. If the limb darkening is taken into ac- 
count in the generated data, the accretion rate obtained with 
the PN model is underestimated by factor ~ 0.8, almost in- 
dependently of the inclination angle. The mass estimation is 
exact at inclination i ~ 45°, underestimated for higher and 
overestimated for lower inclination angles. We stress that 
these fits were performed for a non-rotating black hole only. 

Zhang, Cui and Chen (1997a) investigated accretion 
onto a rotating black hole in an approximate way. They used 
the MCD model and applied GR corrections due to frac- 
tional change in the colour temperature and due to change 
in the integrated flux. As a result, they were able to con- 
strain the black hole spin of a few BHB, including GRO 
J1655-40, which we will analyze in details in Section 3.2 



For comparison, we have applied their method to the spec- 
tra created by our GR model. The result depends on the 
GR corrections applied. Since Zhang et al. (1997) tabulated 
these corrections for a = and 0.998 only, the discrepancies 
are largest for a ~ 0.5. In particular, for spectra generated 
by the GR model with a — 0.998, 0.75 and 0.5, the approxi- 
mate treatment yields a = 0.95, 0.58 and 0.24, respectively. 
We however note, that with GR corrections calculated pre- 
cisely for all spin and inclination angle values the approxi- 
mate method can provide with reasonably accurate results. 

The relativistic disc model parameters are: the black 
hole mass, M, and its spin, a, the inner and outer disc radii. 
Tin and Tout, respectively (in units oi Rg), the accretion rate, 
Md and the inclination angle, i. Subscript 'd' denotes the 
accretion rate in the disc only, excluding all possible energy 
dissipation outside the disc. In our fits we always keep rout ~ 
10'*, so the model has at most 5 free parameters. 

We underline that the overall model spectral shape does 
not change significantly, when these parameters vary. At first 
approximation, there are only two quantities fixing the spec- 
trum: the energy of the peak in the spectrum (or the tem- 
perature) and the total flux (normalization). They translate 
into five model parameters, therefore only two out of them 
can be established uniquely. In most of our fits we fix i, a 
and rin, leaving AI and Md as free parameters. As we show 




Figure 3. Fits of tiie pseudo-Newtonian disc model to the fake 
spectra generated from the GR model. The data were created 
for a non-rotating IOMq black hole with Md = 10^* g s~^ with 
(filled squares) and without (empty squares) limb darkening. 



later, the inner disc radius as a free parameter can be also 
constrained to some extent. 

However, the model, though unable to establish all pa- 
rameters simultaneously, yields tight correlations between 
them. In particular, for a given spectrum and rin ~ rms, as 
a increases, Md decreases (to compensate the increased ac- 
cretion efficiency) but M increases (to keep the absolute in- 
ner disc radius Rin oc VinM when Tin decreases). For a given 
spectrum and a fixed, as rin increases. Aid increases but M 
decreases. Having any additional restrictions for black hole 
mass and inclination angle we can hope to estimate the black 
hole spin. 



3 APPLICATION TO THE DATA 

For spectral fits, we use xspec vlO (Arnaud 1996). The con- 
fidence range of each model parameter is given for a 90 per 
cent confidence interval, i.e., Ax^ ~ 2.7 (e.g. Press et al. 
1992). 

The X-ray spectra of both sources can be decomposed 
into two components: a soft, thermal emission peaking 
around a few keV and a high-energy tail. We interpret them 
in terms of an optically thick cold accretion disc and opti- 
cally thin hot plasma. Therefore we fit the data by a model 
consisting of the GR disc and a high-energy tail, which we 
model either by a power law or by thermal Comptonization 
(see below). The model spectra are absorbed by the inter- 
stellar medium with a column density A^h , for which we use 
the opacities of Balucinska-Church & McCammon (1992) 
and the abundances of Anders & Grevesse (1989). 

While constructing the disc model, we have tabulated 
the transfer function, T, for five values of the black hole 
spin: a = 0, 0.25, 0.5, 0.75 and 0.998. Since the radius of 
the marginally stable orbit, which is the lower limit of rin. 
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is a function of a, we cannot interpolate between transfer 
functions tabulated for different values of a, when rin is close 
to rms- Therefore, we fit the data only for these five fixed 
values of a. In all the fits we keep the spectral hardening 
factor, /col = 1.7 and the limb darkening factor 5 = 2.06. 

3.1 LMC X-1 

3.1.1 Source properties 

LMC X-1 is a luminous X-ray source in the Large Magellanic 
Cloud. It had been a long-standing controversy about the 
optical counterpart of the X-ray source, since the position 
of LMC X-1 established by Emstem was between two stars 
separated by only 6". Cowley et al. (1995) improved the 
position of LMC X-1 from ROSAT-URl observations and 
confirmed that the optical counterpart is a peculiar 07 III 
Star #32 with visual magnitude V ~ 14.8. Using optical 
spectroscopy Hutchings et al. (1987) showed that Star #32 
is in a binary system with orbital period of 4.2288 days. 

They found the mass function, /m = 0.14 ± O.OSM©, 
and the lower limit for the mass of the compact object, 
M > 4M0, which makes LMC X-1 a good candidate for the 
stellar-mass black hole. The inclination angle of the binary 
is constrained by the lack of X-ray eclipses to be i < 63°. 5. 
The mass ratio, q = M,/M is greater than 2, which, together 
with the mass function yields 

Msin^ j > O.8M0. (10) 

The upper limit for the companion star is 25M0, therefore 
M < 12.5M0 and i > 24°. 

The distance to the Large Magellanic Cloud (LMC) is 
uncertain. It is usually published in terms of the distance 
modulus, m—A4 — 5 log Dpc — 5, where m and A4 are appar- 
ent and absolute visual magnitudes and Dpc is the distance 
expressed in parsecs. Different determinations of the dis- 
tance modulus result in conflicting results, from 18.065±0.12 
(Stanek, Zaritsky & Harris 1998) to 18.7 ± 0.1 (Feast & 
Catchpole 1997). This corresponds to the span in distance 
from 38.8 kpc to 57.5 kpc. Since the compact object mass 
determined from the disc observations depends linearly on 
D, this will signiflcantly increase the uncertainty of mass 
determination. 

The X-ray spectrum of LMC X-1 resembles the soft 
state spectrum of Cyg X-1 (Gierlihski et al. 1999). At first 
approximation it can be represented by a blackbody with 
temperature ~ 0.7 keV and a high-energy power-law tail 
beyond the blackbody (e.g. Ebisawa, Mitsuda & Inoue 1989; 
Schlegel et al. 1994; Treves at al. 2000; Schmidtke, Ponder & 
Cowley 1999; Nowak et al. 2001). No transition to the hard 
state has been ever observed. 

3.1.2 Spectral fits 

LMC X-1 was observed by ASCA on April 2-3, 1995. We 
have extracted the SIS spectrum with the net exposure time 
of 23.6 ks. In our fits we use the 0.7-10 keV spectrum, ex- 
cept the channels between 1.8-2.2 keV, which are strongly 
affected by the instrumental gold feature. Since we are in- 
terested in the shape of the disc continuum, this small gap 
in energy channels will not affect our results. 

Both the distance to LMC X-1, D, and the inclination 




E (keV) 

Figure 4. The ASCA/SIS observation of LMC X-1 and a model 
with a non-rotating black hole (upper panel; see Table |). The 
model (solid curve) is decomposed into the GR disc (dotted curve) 
and thermal Comptonization (dashed curve). The lower panel 
shows the data-model residuals divided by errors. 



angle of the system, i, are uncertain. In most of our fits we 
use D = 50 kpc and i = 50°, but we also check our results 
within the acceptable range of 38.8 kpc < D < 57.5 kpc 
and 24° < i < 63°. 5. We do not impose any constraints on 
the black hole mass, M, during the fits, but check its consis- 
tency with the limits obtained from the optical observations, 
4.OM0 < M < I2.5M0, a posteriori. 

For the high-energy tail beyond the disc spectrum we 
chose a thermal Comptonization model (Zdziarski, Johnson 
& Magdziarz 1996). This model is fast and has few param- 
eters, which makes it easy to use. We would like to stress 
that one important difference between the Comptonization 
model and a power law is the low-energy cutoff. The cut- 
off should be present in the Comptonized spectrum around 
the maximum disc temperature if the seed photons come 
from the disc. The soft power law, extending down towards 
lower energies without cutoff, may signiflcantly affect mea- 
surements of the disc parameters. 

The Comptonization model has the following parame- 
ters: the electron temperature, kTe, the asymptotic power- 
law photon spectral index, F, and the seed photons tem- 
perature, kTs. We have found that the electron temperature 
has a negligible effect on the flt results, so we keep it at 50 
keV during our fits. We note, that due to lack of the data 
above 10 keV we cannot constrain plasma parameters, and 
we treat the Comptonization model only as a phenomeno- 
logical component supplementary to the disc model. It is 
not our intention to study the hot plasma properties in this 
paper. 

We flt the data by a model consisting of the GR disc, 
Comptonization and interstellar absorption. First, we as- 
sume that the disc extends down to the marginally stable 
orbit, Tin = rms. The flt results are presented in Table |l[ 
We notice that goodness of the fit does not change with the 
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Table 1. Fit results of the GR disc plus thermal Comptonization model to LMC X-1 ASCA/SIS data for five different values of the black 
hole spin. The inclination angle is i = 50° and the distance to the s ource D = 50 kpc. The disc extends from rjn = rms to rout = 10^ 



The disc and Comptonization models are described in Sections | 
power dissipated in the hot phase). 



and 



3.1.2| , respectively. denotes accretion rate in the disc (excluding 
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Figure 5. Black hole mass versus cosine inclination angle for 
LMC X-1. Curves represent 90 per cent confidence contours for 
fits of the GR disc plus thermal Comptonization model. Distance 
of D = 50 kpc has been assumed. Numbers on contour lines de- 
note black hole spin. The shaded area corresponds to mass and 
inclin ation angle constraints from optical observations (Section 
3.1. l|). 



black hole spin. However, the best fitting mass varies be- 
tween ~ IOMq for the non-rotating black hole to ~ 4OM0 
for the maximally rotating black hole. As we see, only the 
fits with a ^ 0.5 are c ons istent with the upper mass limit, 
M < 12.5M0 (Section |3l|). 

In order to better estimate the black hole spin we fit 
the data in a wide range of inclination angles and compute 
90 per cent confidence contours for inclination-mass rela- 
tion. The results are presented on Figure ^. We find the 
fits with a = 0, 0.25, 0.5 and 0.75 consistent with limits on 
the mass and the inclination angle (shaded area), and only 
a — 0.998 contour lies outside the acceptable area. Even if 
we take into account the significant uncertainty in the LMC 
X-1 distance estimation, this contour is still outside the lim- 
its. Thus, though we cannot precisely find the upper limit 
for the black hole spin, we conclude that the spin close to 
maximal is ruled-out. We stress, however, that this holds 
only when the disc extends down to the marginally stable 
orbit. 

If we allow for free rin > rms we cannot limit the black 



hole spin. In case of a non-rotating black hole the spectral fits 
do not constrain the inner disc radius. The best-fitting value 
is Tin ~ 100, but with the fit improvement by Ax^ = 1-9 only 
in compare to the disc extending down to Vms- This value of 
Tin, however, would require a very small mass of the compact 
object, M ~ 1.5M0 and a super-Eddington accretion rate, 
Md ~ 80LEddC~^. From the lower mass limit, M > AMq, 
we find the upper limit on the inner disc radius, rin ^ 30. 
Since at this radius the black hole spin is not relevant, this 
limit holds for any spin, and we find no constraints on a this 
time. 

Within the acceptable range of M, a, D and i the best- 
fitting relative accretion rate in the disc, rhd = MdC^/^Edd, 
lies between ~ 1 and ~ 5. In terms of the radiated power it 
corresponds to a fraction rjThd ~ 0.06-0.3 of the Eddington 
luminosity, where rj is the accretion efficiency at a given spin. 
The total accretion rate, including the hard tail luminosity, 
is rjTh ~ 0.1-0.5. The unabsorbed bolometric luminosity of 
LMC X-1, inferred from the model, is (2.5-3.7) x 10^® erg 
s"^ (3.2 X 10^* erg s"^ at D = 50 kpc), from which ~ 0.6 is 
in the disc and the rest in the tail. 

The parameters of the Comptonizing tail are worth a 
word of comment. As it can be seen from Table |l], a pho- 
ton spectral index we found, F ~ 3.3-3.5, appears to be 
significantly softer, than typical values of 2.1-2.4 reported 
by Ebisawa et al. (1989) and Schlegel et al. (1994). How- 
ever, the discrepancy emerges from the different model ap- 
plied. When we fit our ASCA data with the simple model 
of the MCD and a power law, we find F = 2.26lo;ig and 
fcTin = 0.807 ± 0.007 keV, which confirms that LMC X-1 
was in the same spectral state during all these observations. 
A different result of softer spectrum was reported by Nowak 
et al. (2001), who fitted the 1996 December 6-8 RXTE ob- 
servation with the multicolour disc + power law model and 
found F ^ 2.9-3.1 and kTin « 0.9 keV. 



3.2 GRO J1655-40 

3.2.1 Source properties 

The X-ray transient GRO J1655-40 was discovered by 
BATSE detector on board CGRO on July 27, 1994 (Zhang 
et al. 1994). The optical counterpart was soon discovered 
by Bailyn et al. (1995a). The inclination of the binary sys- 
tem is in the range 63°. 7 to 70°. 7 (van der Hooft et al. 
1998). The distance to the source, D = 3.2 ± 0.2 kpc was 
derived from jet kinematics (Hjellming & Rupen 1995). The 
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mass function was first obtained by Bailyn et al. (i995b) 
and tlien improved by Orosz & Bailyn (f997), who found 
/iw = 3.24 ± O.O9M0 and classified the companion star as 
F3 IV-F6 IV. In the same paper Orosz & Bailyn determined 
the mass of the compact object with unprecedented accu- 
racy, finding Mx = 7.02 ±O.22M0. However, Shahbaz et 
al. (1999) pointed out that in calculating the radial veloc- 
ity semi-amplitude, Orosz & Bailyn (1997) used observations 
both during X-ray quiescence and during an X-ray outburst, 
which could lead to an incorrect result. Shahbaz et al. (1999) 
used only the X-ray quiescence data, and found different 
value of the mass function, /a/ = 2.73 ± O.O9M0, the mass 
ratio, q = 0.337-0.436, and the compact object mass, M = 
5.5-7.9M0 (95 per cent confidence). With this mass, GRO 
J1655-40 is one of the most firmly established black hole 
candidates. 

Several authors have tried to estimate the black hole 
spin in GRO J1655-40, however they came up with different 
and conflicting conclusions. Zhang et al. (1997a) fitted the 
August 1995 ASCA data of the source with the MOD model 
with relativistic corrections. They found the inner disc ra- 
dius of 23 km for a 7M0 black hole, significantly smaller than 
the marginally stable orbit radius in the Schwarzschild met- 
ric, Rms = 62 km. Their conclusion is that GRO J1655-40 
contains a Kerr black hole rotating at 0.7-1.0 of the maxi- 
mum rate. Sobczak et al. (1999) applied a similar approach 
to soft state RXTE observations, finding average inner disc 
radius of 7?in = 'i.2Rg. They associated this radius with the 
marginally stable orbit corresponding to the spin of a = 0.5. 
Taking into account uncertainties in the spectral hardening 
factor, /col, and the distance to the source, they concluded 
that a < 0.7. Makishima et al. (2000) applied the relativis- 
tic accretion disc model in Schwarzschild metric to the same 
data and found that if the black hole is not rotating, its 
mass, M = 2.9 ± O.IM0, is incompatible with optical con- 
straints, therefore there must be a rotating black hole in 
GRO J1655-40. 

An alternative approach to estimating the black hole 
spin is analysis of quasi-periodic oscillations (QPOs) ob- 
served in the power spectra. Remillard et al. (1999) anal- 
ysed RXTE observations of GRO J1655-40 and found four 
characteristic QPOs. Three of them occupy relatively stable 
frequencies of about 0.1, 9 and 300 Hz. Cui, Zhang & Chen 
(1997) associated the highest frequency QPO with the nodal 
precession frequency of the tilted disc and derived a = 0.95. 
On the other hand, SteUa, Vietri & Morsink (1999) inter- 
preted the highest frequency QPO in terms of the perias- 
tron precession frequency, and the ~ 9 Hz QPO as a second 
harmonic of the nodal precession frequency. This interpre- 
tation yielded a ~ 0.1. Gruzinov (1999) linked the ~ 300 Hz 
QPO with the emission of the bright spot at the radius of 
the maximal proper radiation flux, and inferred a ~ 0.6. 

In X-rays GRO J1655-40 is highly variable and under- 
goes spectral transitions similar to that of classical X-ray 
novae. Ueda et al. (1998) analysed the four ASCA obser- 
vations between August 1994 and March 1996. They dis- 
tinguished four distinct states, named "high", "low", "dip" 
and "off'. The high state was observed by BATSE during 
the outburst in the energy range of 20-100 keV (Tavani et al. 
1996; Zhang et al. 1997b), while the low state corresponds to 
times when the source was weak also in the BATSE range. 
This high state can be associated with the soft state, since 



its spectrum is dominated by the ultra-soft disc component 
with additional high-energy power law. Zhang et al. (1997b) 
found the power- law photon spectral index of F = 2.43 ±0.3 
from the BATSE data simultaneous to the ASCA observa- 
tion analysed in this paper. A long-time RXTE monitoring 
(Sobczak et al. 1999) shows that the soft state can be de- 
scribed by the MCD model with fcTin ~ 0.7 keV plus a power 
law with photon spectral index of F ~ 2-3. 

3.2.2 Spectral fits 

GRO J1655-40 was observed by ASCA on August 15-16, 
1995. This is the observation from epoch HI from Ueda et al. 
(1998), who associated it the high X-ray state of the source. 
From this observation we have extracted the GIS spectrum 
with the net exposure of 3810 s after dead-time correction. 
We use the data in the 1.0-10 keV range. 

We assume the distance to the source, D = 3.2 kpc 
(Hjellming & Rupen 1995), and the inclination angle of 
70° (Orosz & Bailyn 1997; van der Hooft et al. 1998), un- 
less stated otherwise. We do not impose any constraints on 
the black hole mass during the fits, but check its consis- 
tency with the limits obtained from the optical observat ions , 
5.5M0 < M < 7.9M0 (Shahbaz et al. 1999; see Section |!|), 
a posteriori. 

The high-energy tail beyond the disc spectrum is less 
significant than in the case of LMC X-1 (see Figures ^ and 
and we find that the particular choice between a Comp- 
tonization model and a power law does not affect the disc 
fit results, so we chose a power law with the photon spectral 
index fixed at F = 2.4 (Zhang et al. 1997b). On the other 
hand, the tail cannot be neglected. A fit without the tail 
is worse by Ax^ = +90 (at 188 d.o.f.) in compare to our 
best fit (see below) and shows increasing positive residuals 
towards higher end of the energy band. 

We fit the data by a model consisting of the GR disc, 
power law and interstellar absorption. First we assume the 
disc extends down to the marginally stable orbit, rin — rms. 
The model fits the data fairly well, with ~ 230/187, how- 
ever there is a strong residual pattern below ~ 3 keV (see 
panel (b) on Figure h). The pattern is likely to be gener- 
ated by atomic processes. However, the interstellar absorp- 
tion, both on neutral and ionized matter, cannot explain 
this pattern. Varying elemental abundances in the absorber 
does not improve the fit. Also, the pattern cannot be fit- 
ted by ionized Compton reflection. We notice that relativis- 
tically broadened edge is required to explain the observed 
spectrum. 

Therefore, we included a smeared absorption edge 
(smedge model in xspeC; Ebisawa 1991) in our model spec- 
trum, smedge is only a phenomenological formula invented 
to reproduce the relativistic smearing and its parameters 
should not be taken as real physical quantities. In particular, 
when the smearing is introduced the model edge threshold 
energy is lower than the rest frame threshold energy of the 
physical edge. We have found that inclusion of the absorp- 
tion edge improves the fit considerably, yielding much bet- 
ter ~ 165/185. The smearing width of the edge is poorly 
constrained around ~ 0.5 keV, so we fix it at this value dur- 
ing all fits. The threshold energy, typically 1.25 ± 0.05 keV, 
might be consistent with the Ne-like Fe L-shell absorption 
edge (at 1.26 keV), originating in the ionized disc. Detailed 
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Figure 6. The ASCA/GIS spectrum of GRO J1655-40 and the 
best-fitting model with a = 0.75 (see Table Panel (a) shows 
the data and the model. The model (solid curve) is decomposed 
into the GR disc (dotted curve) and the power law (dashed curve). 
The residuals (data-model divided by errors) shown in panel (c) 
correspond to this fit. For comparison, on panel (fe) we present 
the residuals of the fit without absorption edge. 



analysis of this feature is however beyond the scope of this 
paper. 

Finally, our best-fitting model consists of the GR disc, 
the power law, the interstellar absorption and the absorption 
edge. Fit results for five values of the black hole spin are 
presented in Table ^. The data and the best-fitting model 
for a — 0.75 are shown on Figure ^ 

As it was in the case of LMC X-1, spectral fits do not 
constrain the black hole spin. However, the spin can be con- 
strained taking into account mass limitations. Only the fit 
with a = 0.75, yielding M = 5.84lo;o5M0, is consistent with 
the mass constraint of (5.5-7.9)Mq , inferred from the opti- 
cal observations by Shahbaz et al. (1999). Again, we fit the 
data in a range of inclination angles and compute 90 per 
cent confidential contours for inclination-mass relation. The 
result is presented on Figure |^. Due to better quality of the 
data and due to the fact that GRO J1655-40 spectrum is 
dominated by the disc emission, these contours are much 
narrower than in the case of LMC X-1. Only the a = 0.75 
contour overlaps with the shaded area indicating allowable 
values of the inclination angle and the mass. More precise 
limits on the spin can be found by interpolating between the 
values found from the fits. If the GR effects are neglected, 
for the given spectrum the disc must keep constant inner ra- 
dius while varying M and a. Providing that rin = 7"ms, this 
leads to a simple relation between the mass and the spin 
M oc r-~s(a), where rms(a) is given by equation The GR 
effects modify this relation and we have found that a func- 
tion M = Ci -\- C2r~s (a) approximates the data well, so we 
use it for interpolation. Taking into account uncertainties 
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Figure 7. Black hole mass versus cosine inclination angle for 
GRO J1655— 40. Curves represent 90 per cent confidence contours 
for fits of the GR disc plus power law model. Distance ol D = 
3.2 kpc has been assumed. Numbers below contour lines denote 
the black hole spin. The shaded area corresponds to mass and 
inclin ation angle constraints from optical observations (Section 
3.2.l|). 



in the distance to the source, we have found that for given 
limits 5.5 < M/Mq < 7.9 and 64° < i < 71° the black hole 
spin is 0.68 < a < 0.88. In this range of spin values return- 
ing radiation is negligible (see Section]^, and does not affect 
our result. 

Next, we fit the data allowing the inner disc radius, rin, 
to be a free parameter. The results are presented in Table ^ 
In all fits we find that the inner disc radius is consistent with 
Tms, however larger values are possible. Now the a — 0.998 
fit became acceptable with M = lOt^ ^M©. Since the upper 
limits for the mass remained virtually the same after freeing 
Tin (see Tables ^ and |^, the lower limit on the spin also 
remained the same, a > 0.68, this time. 

Within the acceptable range of M, D and i (assuming 
= Tyna) the relative accretion rate, rhd, is between ~ 0.6 
and ~ 1.9. In terms of the radiated power it corresponds 
to a fraction r^rn^ ~ 0.1-0.2 of the Eddington luminosity. 
The unabsorbed bolometric luminosity of the disc in GRO 
J1655-40, computed from the model, is (1.2-1.6) x 10^** erg 
s"^ (1.4 X 10^* erg s"^ at D = 3.2 kpc). In the observed en- 
ergy range only a small fraction of the luminosity is radiated 
in the tail. 



4 DISCUSSION AND CONCLUSIONS 
4.1 Black hole mass and spin 

Optical and UV observations of black hole binaries can pro- 
vide with the mass function, which in turn can help estimat- 
ing, to some extent, the black hole mass and the inclination 
angle of the binary. We know about a dozen of X-ray Galac- 
tic binaries for which the mass function implies the compact 
object mass larger than 3M0, making them good candidates 
for black holes. However, the black hole spin cannot be es- 
timated this way. 
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Table 2. Fit results of the GR disc plus power law and an absorption edge model to GRO J1655-40 ASCA/GIS data for five different 
values of the black hole spin. The inclination angle is i = 70° and the distance to the source D = 3.2 kpc. The power-law photon index 
is fixed at r = 2.4. The disc extends from rin = Vms to rout = Iff*. The disc model is described in Section jij 



a 


0.00 


0.25 


0.50 


0.75 


0.998 


Nil (10^^ cm-2) 
Sedge (keV) 

''"edge 

M (Mo) 

Md (IQis g s-i) 


'■^"-0.18 

1 9t:+0.05 

n 1O+0.03 

2 65+0 03 
^•D3_0.02 

O r,i+0.04 

■'■^^-0.03 


' -^"-0.17 
, 9C+0.05 
-^■^^-0.05 
13+004 

9 71 +0.04 
^- ' -^-0.04 


7 50+0-17 
^ 9C+0.05 
^■^^-0.05 
14+0-04 
^■^^-0.03 
4 07+0-04 

9 1C-+0.03 
^•^^-0.02 


'•='^-0.18 
1 9r-+0.04 

16+"-°'' 
'^•^"-o.os 

^ O.+0.06 

1 CI +0.02 
^■^^-0.02 


7 70+0-15 

' ■ "-'-0.17 
I 90+0.04 
■ '^-0.04 

]^7+o.03 

-0.03 
16 2+°-l 
560+0-007 


X^/ISS d.o.f. 


166.1 


165.3 


164.3 


163.4 


165.0 



Table 3. The same as in Table ul, but now rin is a free fit parameter. For the sake of clarity, we show only rin and M. 
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In this paper we have shown an application of the GR 
disc model to the soft state X-ray spectra of LMC X-1 and 
GRO J1655-40 and demonstrated how it can restrict the 
black hole spin. The model depends on five basic parameters: 
the black hole mass and spin, the disc inclination angle, 
the accretion rate and the inner disc radius. Generally, the 
spectral shape changes only slightly with these parameters, 
and only two out of them can be established uniquely by 
fitting the observational data. Additional assumptions about 
some of the parameters are necessary. In most of our fits we 
have fixed i and assumed rin — rms- Then, for several fixed 
values of a we have successfully fitted M and M. Results 
presented in Tables ^ and |2| show that remained virtually 
the same through the whole range of a and the black hole 
spin cannot be reckoned this way. On the other hand, at a 
given spectral shape there is strong correlation between the 
parameters, in particular between M and a when i and rin 
are fixed. With an independent estimation of the mass and 
the inclination angle, we can constrain the black hole spin. 
The tighter constraints on the mass and the inclination, the 
more accurate spin estimation can be obtained. 

The mass and the inclination angle of LMC X-1 are 
evaluated rather poorly. Therefore, it is not possible to give 
a good constrain on the black hole spin. We can only rule out 
a black hole which is close to, or maximally rotating. GRO 
J1655-40 gave us better chance. The mass and the inclina- 
tion angle are measured with relatively high accuracy. Also, 
due to better quality of the X-ray data, we have obtained 
smaller statistical errors on M. As a result, we could limit 
the spin, 0.68 < a < 0.88 (when rin = rms). This result is 
only weakly affected by neglecting returning photons, which 
are negligible for spin less then 0.9. 

We underline, that when we allow for free rin > rms, no 
upper limit for the black hole spin can be imposed. This is 
because for a given spectral shape increase of Hn with fixed 
M and i leads to increase of a. However, there are clues that 
the accretion disc in BHB extends down to the marginally 
stable orbit all the time in the soft (high or very high) state. 

Long-term 1996-1997 RXTE monitoring of GRO 



J1655-40 (Sobczak et al. 1999) shows remarkable constancy 
of the MOD inner disc radius over the period of a few hun- 
dred days, while the source went through the high and very 
high spectral states. There is an exception of a few ob- 
servations during the very high state, where the inner ra- 
dius dropped suddenly by factor ~ 3. Merloni, Fabian & 
Ross (2000) interpret it as a drop in the disc accretion rate, 
which caused increase of the hardening factor and decrease 
of the apparent inner radius, while the actual rin has not 
changed. If the inner disc radius remained constant indeed, 
it is reasonable to assume that the disc extended down to 
the marginally stable orbit all the time when the source was 
in the high or very high state, including also the ASCA ob- 
servation analysed in this paper. 

Long-term 1996-1998 RXTE monitoring of LMC X-1 
(Wilms et al. 2001) shows in turn significant variations of 
the MCD inner disc radius. However, it is not clear whether 
these variations can be attributed to real changes of rin. 
They might be, at least in part, due to systematic errors in 
the model. 

So, does Tin vary during the soft state or not? If the 
cold disc in either of the above observations was indeed trun- 
cated, it should turn into an optically thin inner fiow below 
Tin. Hot plasma in the inner flow could be a source of the 
high-energy tail beyond the disc spectrum. However, opti- 
cally thin solution can exist only for luminosities lower than 
a few per cent of the Eddington luminosity. At higher accre- 
tion rates (and higher densities) Coulomb transfer of energy 
from the protons to the electrons becomes efficient, the op- 
tically thin flow becomes radiatively efficient and collapses 
to the optically thick Shakura-Sunyaev disc (Chen et al. 
1995; Esin, McClintock & Narayan 1997). In both observa- 
tions analysed in this paper luminosity exceeded 10 per cent 
of Lfidd, therefore we should not expect any hot, optically 
thin flow. Instead, we propose that the cold disc extends all 
the way down to the marginally stable orbit, and the high- 
energy tail emission is produced in active regions above the 
disc. Hence, we find the assumption of rin = rms as well 
founded, and our best constraint on the black hole spin in 
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GRO J1655-40 is 0.68 < a < 0.88 and in LMC X-1 we rule 
a black hole rotating with the spin close to maximum. 

The most spectacular feature of GRO J1655-40 is its 
radio jets (Tingay et al. 1995; Hjellming & Rupen 1995), ob- 
served also in another high-spin black hole candidate, GRS 
1915+105 (Mirabel & Rodriguez 1994). Co-existence of rel- 
ativistic jets and rotating black holes in BHB makes and in- 
teresting link to the problem of radio dichotomy of quasars, 
and supports the black hole "spin paradigm", according to 
which jets in quasars are powered by rotating black holes 
(e.g. Moderski, Sikora & Lasota 1998). Since no strong ra- 
dio jets were found in LMC X-1, the non-rotating black hole 
should be expected in this system. Unfortunately, our re- 
sults do not allow us to make such claims. Further system- 
atic studies of several jet and non-jet sources could provide 
with more evidences for the spin paradigm. 



4.2 Compton reflection 

The X-ray spectra of LMC X-1 and GRO J1655-40 are sim- 
ilar to the soft state spectrum of Cyg X-1 (Gierliiiski et al. 
1999), though the first two sources are significantly brighter. 
An important difference is that we do not find Compton re- 
flection features here, while Cyg X-1 in the soft state showed 
strong reflection component with covering angle ^/2n ~ 0.7 
and a broad Fe Ka line. This can be explained by the in- 
strumental limitations or/and by the intrinsic difference in 
the source geometry and energetics. The Compton reflection 
from the cold matter can be identified by Fe Ka features 
around 7 keV and by reflection continuum, which peaks 
around 30 keV in the EFe spectrum. Since ASCA can ob- 
serve only up to ~ 10 keV, the detection of the reflection 
continuum with ASCA data only might be difficult, if not 
impossible. The iron K edge and line parameters are sensi- 
tive to the shape of the underlying continuum, which can- 
not be properly established without the high-energy data. 
For example, Dotani et al. (1997) analysed the ASCA data 
of Cyg X-1 and did not find any Fe line. Gierliiiski et al. 
(1999) re-analysed the same data jointly with the simulta- 
neous RXTE observation, and found the presence of the 
iron line with high statistical significance. LMC X-1 spec- 
trum has poor statistics above ~ 5 keV, so the reflection 
features can escape the detection. The high-energy tail in 
GRO J1655-40 is much weaker than in Cyg X-1. At 7 keV 
the disc radiates ~ 3.5 times more energy than the tail. 
Therefore, the Fe K-shell features will be very weak in the 
total spectrum, and difficult to detect. In the soft state of 
Cyg X-1 the situation is opposite: the disc emission at 7 keV 
is negligible (see Figure 8 in Gierliiiski et al. 1999 and Figure 
^ in this paper). Moreover, the Fe K-shell features coming 
from the disc around the fast-spinning black hole are signif- 
icantly smeared, which additionally hinders the detection. 

We note, that ASCA and RXTE detected characteristic 
iron features in GRO J 1655-40, at different periods (Ueda et 
al. 1998; Tomsick et al. 1999; Balucihska-Church & Church 
2000; Yamaoka et al. 2000). However, precise ASCA/ RXTE 
observations show that these features can be consistently ex- 
plained by resonant absorption lines in the corona above the 
disc without reflection components (Yamaoka et al. 2000). 



4.3 Disc stability 

An unclear issue is stability of the cold disc. Both sources 
are very bright with luminosities being a significant frac- 
tion of the Eddington luminosity (rhd ~ 1-5 is a very con- 
servative limit for LMC X-1). A radiation pressure dom- 
inated region, which is thought to be unstable, arises in 
the disc above the critical accretion rate, merit. In the 
pseudo-Newtonian approximation, for a lOM© black hole 
and a — 0.1, rhcrit ~ 0.64 (Gierlinski et al. 1999; taking into 
account accretion in the cold disc only, i.e. / = 0). Therefore, 
the discs in both sources seem to be dominated by the radia- 
tion pressure. Since we consider here rotating black holes, we 
checked this result by calculating the ratio of the radiation 
pressure to the gas pressure in the Kerr metric (Novikov & 
Thorne 1972), for all acceptable fit parameters in Tables ^ 
and |2| We found that in both sources the inner part of the 
cold accretion disc below several tens of Rg is dominated by 
the radiation pressure indeed. Thus, both discs should be 
unstable against secular and thermal instabilities. However, 
despite this inconvenience the discs in LMC X-1 and GRO 
J1655-40 apparently do exist. 

This riddle probably arises from our lack of understand- 
ing of microscopic mechanisms of accretion. The instability 
develops in the standard a-prescription theory, where r(j>- 
element of the viscous stress tensor, Tr,p, is proportional to 
the total pressure, i.e. to the sum of the gas pressure and 
the radiation pressure. The viscosity parameter a includes 
not well understood physics of viscosity, which is supposed 
to be due to chaotic magnetic fields and turbulence in the 
gas flow. It is not quite clear why T^^ should be proportional 
to the total pressure. If we assume that T^^ is proportional 
to the gas pressure only (the so-called /3-disc; see e.g. Stella 
& Rosner 1984), the disc becomes stable in the radiation 
pressure dominated region (Figure H). The existence of the 
radiation pressure dominated discs should impose significant 
constraints on the theoretical models of viscosity in the ac- 
cretion flow. 



4.4 Model reliability 

One important question we should ask is reliability of the 
model and parameters inferred. While modelling the disc 
spectrum in the "proper" way one should solve the following 
four problems, (i) The radial disc structure, (ii) the vertical 
disc structure, (iii) the radiation transfer and (iv) the rela- 
tivistic effects. Our understanding of the first three problems 
is poor and there are still a lot of open questions concern- 
ing the viscosity prescription, vertical distribution of gravi- 
tational energy dissipation, effects of radiation pressure and 
irradiation of the disc, importance of bound-free opacity and 
heavy elements, and importance of Comptonization, just to 
name the most important issues. Consequently, every single 
prediction available in the literature is based on simplifying 
assumptions and is model dependent. On the other hand, the 
simplest solution, i.e. the sum of blackbodies gives the best 
description of the observational data. For example, a smooth 
multicolour spectrum with no ionization edges consistently 
reproduces the AGN observations in which the Lyman edge 
(expected by more advanced models) is not observed. As we 
have shown in this chapter, the multicolour model fits the 
black hole candidate spectra perfectly as well. 
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Figure 8. Two solutions of the optically thick disc around a 
IOMq non-rotating black hole at the radius of the most significant 
instability of r = 15.1. The viscosity parameter a is 0.1. The left- 
hand curve represents the solution with T^.^ ex Pjad + ^gas (a- 
disc). The right-hand curve shows the solution with T^^ oc Pgas 
(/3-disc). The shaded region shows the most conservative limits 
on the disc accretion rate of LMC X-1 at a = 0. The dashed 
line represents the Eddington luminosity (efficiency rj = 1/16). 
The inner disc is dominated by the radiation pressure, where the 
o-disc solution is unstable and the /3-disc solution is stable. 



We correct the emerging spectra for two effects that 
significantly affect measurements of the disc normahzation 
and temperature. The first one is a deep gravitational po- 
tential in the vicinity of a fast rotating black hole. All special 
and general relativistic effects are treated to high accuracy 
by means of the transfer function, as was described in Sec- 
tion ^. The second important aspect is the effect of electron 
scattering, which we approximate by a single colour tem- 
perature correction factor, /coi, using a diluted blackbody 
as a local spectrum (see Section P). Under certain simplify- 
ing assumptions Shimura & Takahara (1995) computed the 
vertical disc structure around a Schwarzschild stellar-mass 
black hole and solved the radiation transfer. They found that 
for the viscosity parameter a = 0.1 the local spectrum can 
be represented by a diluted blackbody when Comptoniza- 
tion is effective, which takes place in the inner disc region 
when md ^ 1. For M = lOM© and rh^ — 1 they found that 
the sum of diluted blackbodies fits the whole disc spectrum 
very well with the best-fitting hardening factor of /coi = 1.7 
(see Figure 2 in Shimura & Takahra 1995). Merloni et al. 
(2000) made a similar comparison between the multicolour 
disc and a more realistic disc model with the vertical tem- 
perature structure and radiation transfer solved and came 
up with a comparable result. The multicolour disc spectra 
fit the realistic disc spectra well, and for md ^ 1 the hard- 
ening factor is /coi ~ 1.8. In both observations in this paper 
rhd ~ 1, therefore, within the trustworthiness limits of the 
Shimura & Takahara and Merloni et al. calculations, we find 
our model reliable. A small variation in the hardening factor 
does not change our results significantly. For example, the 
GRO J1655-40 spin for /coi = 1.8 is 0.60 < a < 0.83. 
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APPENDIX A: CONSTRUCTION OF THE 
PHOTON TRANSFER FUNCTION 

We consider a gravitational field of a rotating black iiole 
with the mass M, and the angular momentum, J. In both 
appendices we use G — c — 1 units, the ( — h ++) signature 
and the Boyer-Lindquist coordinates, {t, R,0, (j)). Non-zero 
components of the metric tensor are given by 



gtt 


= -(l-2r/S), 


9t4. 


= -2arsin^6'/S, 




( 2 ^ 2 ^ 2a- 


g<t><t> 


= 1 r -f a H 




= E/A, 


gee 


= E, 



where 



A = + - 2r, 
Yi — + cos^ 9 



(Al) 



(A2) 



and we use dimensionless distance and specific angular mo- 
mentum parameter 

Our approach to calculation of the disc photon trajec- 
tories mostly follows that of Cunningham (1975). A null 
geodesic in the Kerr metric is described by two constants of 
motion, A and 5, defined as 

L .2 g 

EM' 



A 



(A4) 



where E is the photon energy at infinity, L is projection of 
the photon angular momentum on the symmetry axis and 
C is the Carter constant (in the definition of ^ we took 
into account the condition C > satisfied by trajectories 
intersecting equatorial plane.) 

We take an usual assumption that the radial velocity 
in the disc can be neglected and its motion is Keplerian. In 
the non-rotating frame, the angular velocity and the linear 
velocity, respectively, are given by 

' ' (A5) 



K 



M r3/2 + , 



and 
Vk = - 
where 



/ 2ar 



A 4 , 2 2 . o 2 

A — r + a r + 2a r, 



(A6) 



(A7) 



and UK = MQk ■ Then, we get the following relation between 
the constants of motion and the emission angles in the disc 
rest frame (e.g., Cunningham & Bardeen 1973) 

sin 8e sin (j}e + Vk 



A = 



C = 



(ri Ai/2 -I- 2are Vk)M + sin de sin < 



(1- Vk)"'''^(1-Ai.^k)cOS Oe 



(A8) 



where is the initial photon radius, 9^ is a polar angle 
between the photon initial direction and the normal to the 
disc and 0e is the azimuthal angle, in the disc plane, with 
respect to the r-direction. 

Then, the redshift of the photon emitted from the disc 
at a distance re is given by 



gcff = r< 



I^Ay/2 (i_l/i^)i/2 



1 X ' (^9) 

1 — (jJkA 

and the angle, i, at which the photon will be observed far 
from the disc in the flat space-time is determined by the 
integral equation of motion 



5R 



'1/2 



dr = 



-1/2 



d9. 



(AlO) 



Tv/2 



where K and 6 are radial and polar effective potentials, 
respectively, 

AVsin^ 9) . 



K(r) 
Q{9) 



(r 



^2 , 2 , 

: ^ + COS I 



(All) 



The photon transfer function is tabulated in discrete 
steps in five dimensions (a, re, cosz, g^g, cos^e). Generated 
photons are summed in the appropriate elements of the pho- 
ton transfer function table according to the following algo- 
rithm. 
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(i) for given a and generate the photon initial direction 
from the distribution uniform both in 6^ and cos 6e , 



(ii) find the constan ts of motion from equation (A8), 

(iii) solve equation (AlO) for i (initial sign of 5R in equa- 
tion (AlO) is negative for n/2 < (ji^ < 37r/2 and positive 



otherwise) 

(iv) find Qcs from equation (A9), 

(v) increment an element of the transfer function table 
corresponding to a, re, cosi, g^a and cos^e- 

In order to compute and tabulate the transfer function we 
have calculated trajectories of 10^ photons. The disc is axi- 
ally symmetric and emission from a given point of the disc 
is also axially symmetric. Therefore, integration over 0e for 
trajectories starting at a given point of the disc and observed 
at a given angle i and collected at all angles cj)o is equivalent 
to integration over all trajectories starting at a given radius 
re, which are observed at a given angle i and any particular 



Ps = e;: ( - e% ef /el 



4'' 



(B7) 



The factor 0$ gives correction for the integral element drd(f> 
in equation (|l|) (integration over cj> is hidden in the convo- 
lution) which should be taken into account in calculation of 
the number of photons emitted from the accretion disc. 



APPENDIX B: TRANSFORMATION OF THE 
PHOTON NUMBER INTENSITY 

As pointed out in Section |^, application of the photon trans- 
fer function requires transformation of the photon number 
intensity between the reference frames of the distant ob- 
server (the Boyer-Lindquist coordinate frame) and the ob- 
server co-rotating with the disc. The latter one is represented 
by the covariant tetrad (we show only the components used 
below): 

X {gtt + <^K5t0, 0, 0, gt4, + c-JKg4,4,) , 
e"-' = (0,r/VA,0,0), 

""^ — r — ^ (a") " i-^ytt ~r yt<p 



VkQu + gt4, 



2arVK\ „ „ /r^A^/^ 2arVK\ 
H — j ,0,0, VKgt,t> + gu I — \ — I 



(Bl) 



where order of vector components is given by {t,r,6,(f>), 
and primes denote coordinates defined in the disc local rest 
frame. Then, the one-form At' can be written in terms of the 
Boyer-Lindquist coordinates as 

At' = EtAt + e^Afj) (B2) 
As for the observer at rest in the disc A(f)' = 0, and 
d<j!>' = ef dt + e^'d<j!>, (B3) 
we obtain the following relation 

At' = PtAt, (B4) 
where 

Pt = el - e^ef /e^ . (B5) 
Similarly, for the disc unit area we find 

dr'd0' = e^'dr (^ef At + ef df/.) = /3sdrd0 (B6) 

(in the latter equality we have used the condition At' = 0, 
relevant for the three-space measurements in the disc rest 
frame), where 
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